Spectral portraits of the Orr— Sommerfeld operator with 

large Reynolds numbers. 
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^ ' Abstract. A model problem of the form 

O-i. 

^ ! -iey" + q{x)y = Xy 

2/(-i) = j/(i) = o 

I is associated with well-known in hydrodynamics Orr-Sommerfeld operator. Here A is the spec- 

tral parameter, e is the small parameter which is proportional to the viscocity of the liquid and 
to the reciprocal of the Reynolds number, and q{x) is the velocity of the stationary flow of the 
' liquid in the channel ^ 1. We study the behaviour of the spectrum of the corresponding 

CO I model operator as e —> with linear, quadratic and monotonous analytic functions. We show 

' that the sets of the accumulation points of the spectra (the limit spectral graphs) of the model 

I and the corresponding Orr-Sommerfeld operators coincide as well as the main terms of the 

' counting eigenvalue functions along the curves of the graphs. 

o ■ 

p ^. The well-known in hydrodynamics Orr-Sommerfeld equation is obtained by linearization of 

the Navier-Stokes equation in the infinite three-dimension spatial layer {x,^,ri) G M.^, where 
|x| ^ 1 and rj) G M^, assuming that a stationary unperturbed solution for the velocity profile 
is of the form {q{x), 0, 0). This equation with respect to a function y = y{x) has the form (see 
details in the book of Drazin and Reid for example) 



> 



(0.1) {D^ ~ayy-iaR[qix)iD^ -a^)~q"ix)]y = -laRXiD^ -a^)y. 

I Here D = d/dx, a is a wave number (a ^ 0) appearing after the separation of the variables 
(^,77) G M^, R is the Reynolds number characterizing the viscosity of the liquid and A is the 
spectral parameter. Usually the boundary conditions 

(0.2) y{±l) = y\±l) = 

are associated with equation (|0.1|) . 

The main goal of this paper is to describe qualitatively the limit behaviour of the spectrum 
of problem (jO-H) . ()0.2|) as i? ^ 00. The Reynolds number R is proportional to the reciprocal 
of the viscosity. Therefore, our problem is devoted to the description of the spectrum of the 
Orr-Sommerfeld problem for the liquid which is almost ideal. During long time it was supposed 
that the spectrum of the Rayleigh problem 

qix){D' - a')y - q"ix)y = \{D' - a')y, 
2/(-l) = y(l) = 0, 
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plays an important role for the solution of the problem which is raised above. The Rayleigh 
problem is obtained (after dividing of equation (|(J.1|) by —laR) by formal passing to the limit 
as i? — s> oo and by eliminating of the "superfluous" boundary conditions. A vast literature 
is devoted to the study of the spectrum of the problem (j(J.3|) . The papers of Lin [B] and the 
cited book [2j can be recommended as a guide in this subject. Actually, the spectrum of 
the Rayleigh problem plays no essential role in the description of the spectral portraits of the 
Orr-Sommerfeld problem as ^ oo. This will be seen in the sequel. 

It is known [5] , that the spectrum of the problem ()0.3|) consists of the segment [m, M] , where 
m and M are minimum and maximum of the function q{x) (it is supposed that q{x) is contin- 
uous), and possibly, some isolated eigenvalues of finite multiplicity lying outside this segment. 
Perhaps, Heisenberg was the first who noticed that there is no continuity between the spectra of 
the Orr-Sommerfeld problem with large R and the spectrum of the Rayleigh problem. It may 
happen that there is a domain containing the interval (m, M) which is free of the spectra of 
problem (jO.H) . ()U.2|) for all sufficiently large R. This phenomenon was named as a "Heisenberg 
tongue". In 1924 Heisenberg proved the existence of a system of fundamental solutions for 
equation ()0.1|) having special representation (see [Ij). This result is very essential for explana- 
tion of such a phenomenon. However, we are not aware of Heisenberg' papers containing ideas 
to explain the phenomenon. It was Morawetz T], who proved that in the case of the Couette 
profile q{x) = x the spectrum of the problem (|(J.1|) . (|(J.2j) is concentrated in a (5-neighbourhoods 
of the ray [—i/ \/3, ~ioo), of two segments [±1, —i/\/3] and the isolated eigenvalues {/ifc} of the 
Rayleigh problem for q{x) = x. Actually, the Rayleigh problem for q{x) = x has no isolated 
eigenvalues (see pi), therefore, the last reservation can be droped. Moreover, Morawetz made 
an attempt to prove a similar result for more general functions. She assumed (this assumption 
is quite essential in her method) that q{x) is an entire function with real values on the real 
axis and maps bijectively the whole complex plane C into itself. Probably, it was not noticed 
that all the function possessing these properties have the representation q{x) = Qx + 6, where 
Q, 6 E M.. There is another important problem which was left beyond the paper |7j. What 
is the set of the accumulation points of the eigenvalues as i? — *• oo? Morawetz emphasized in [7] 
that her method did not allow to get the information if the eigenvalues exists as i? ^ oo near 
each point of the segments [±1, 

In 90-ies there appeared papers (see \T7\ . |10j . for example), where a simplier problem of 
the form 

(0.4) -iez" + q{x)z = Xz, 



was associated with (jO.H) . ()0.2|) . Here e and A are the small and spectral parameters, respec- 
tively. One can consider this problem as a simplified model for (jO.H) . ()0.2|) . The following 
arguments can confirm this point. 

Let us make in equation ()0.1|) the substitution z = {D"^ — a'^)y. From this equation taking 
into account the conditions y{—l) =?/'(— 1) = 0, we find 



(0.5) 



z{-l) 



X 



(0.6) 




-1 



Then equation (j0.1|) takes the form 
(0.7) -ie{D^ 



c?)z + q{x)z + Kz = \z, 
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where 



X 

Kz = j sm\i2{x-i)q"{i)z{i)di, 



1 

aR 



-1 



It follows from that y(— 1) = y'{—l) = 0. Therefore, one can rewrite boundary condi- 
tions ()().2j) in the form 

1 1 
(0.8) /.«)si„ha(l-O^^0. /.K)cosha.(l-O<i^^0. 

-1 -1 

Hence, problem (jU.ip . (j(J.2|) is equivalent to problem (|U.7p . (j(J.8|) . This reduction was realized by 
Orr in 1915. Now, if we neglect the influence of the integral operator K (notice, that = in 
the case q{x) = x) and assume that boundary conditions do not change essentially the spectral 
portraits as e ^ 0, we come to model problem ()0.4p . (jO.Sp up to the shift of the spectral 
parameter by iea"^. Of course, these arguments are heuristic. However, we shall show below 
that the similarity of the spectral portraits of the model and the Orr-Sommerfeld problem as 
e — >■ can be proved rigorously. 

If we put £ > in equation ()0.4p instead of ie, we get a selfadjoint problem with small 
parameter. Such a problem is well investigated long ago (see |12j . for exemple). Its spectrum 
is real, it condenses as e ^ and the formulae for the eigenvalue localization can be written 
out explicitly. They are known as the Bohr-Sommerfeld quantization formulae. 

The replacement of e by ie changes the problem dramatically. In 1997 the author |13j 
described the spectral portraits of model problem ()U.4|) . ()0.5|) as £ —> in the case q{x) = x 
and paid attention to the fact that for analytic functions q{x) the spectrum is concentrated 
along some curves which are determined by the geometry of the Stokes lines of equation ()0.4|) . 
It is useful here to recall the definitions. The zeros of the equation g (2;) — A = in the complex 
2;-plane are called turning points and the lines 



= i 2; G 



Re / Vz(g(O-A)rfe = 



outgoing from a fixed turning point are called the Stokes lines of equation ()().4|1 . These lines 
come either to the boundary of a domain G, where q{z) is holomorphic (in particular, go to 
00 if q{z) is an entire function), or end in other turning points. The maximal connected set 
consisting of the Stokes lines and containing a turning point is called the Stokes complex 
corresponding to The Stokes complex may contain several turning points but not necessary 
all of them. The union of Stokes complexes corresponding to all turning points is called the 
Stokes graph. 

As we have mentioned before, for analytic (or piece- wise analytic) functions q{x) the spec- 
trum of problem ()().4|1 . ()().5|1 is concentrated near some curves in the A-plane as e ^ 0. We call 
them limit spectral curves. The union of the limit spectral curves we call limit spectral graph. 
One should not mix the limit spectral curves or the limit spectral graph with the Stokes lines 
and the Stokes graph which lie in z-plane. Concrete forms of the limit spectral curves can be 
found comparatively easy only for the linear function q{x) = x. For nonlinear functions q{x) 
these curves take a complicated form. Therefore, attempting to describe the spectral portraits 
for non- selfadjoint problems with small or large parameter one meets serious difficulties even 
in the case of model problem ()0.4|) . ()0.5|) . Now we can formulate our goals more concretely. 



SPECTRAL PORTRAITS 



4 



(1) To find functions q{x) of particular and general form for which the spectral portraits of 
problem (j(J.4j) . (jU.Sp as e ^ can be described completely. Of course, special attention 
should be paid to profiles which correspond to stationary solutions of the Navier-Stokes 
equation, in particular to the Couette profile q{x) = x, Poiseuille profile q{x) = 1 — 
and Couette-Poiseuille profile q{x) = ax^ + bx + c, a,b,c & M; 

(2) To find the formulae for the eigenvalues near the limit spectral curves as e 0, 
provided that these curves are already determined; 

(3) To solve the same problems for the original Orr-Sommerfeld problem with the same 
functions q{x). 

At the present moment the problems, which we have raised up, are solved only partially. In 
this paper we present the obtained results based on the research of the author and his students 
A. V. Dyachenko, S. N. Tumanov and M. I.Neiman-zade. We should also mention recent papers 
of Chapman |T], Redparth and Stepin |16j which are closely connected with the problems 
in question. 

1. Model problem: the case q{x) = x. 
Let us consider the spectral problem 

(1.1) -^ey" = {x-X)y, 

(1.2) y{-l) = y(l) = 0, 

where £ is a small parameter. As we have mentioned, computer calculations give a striking 
result: the eigenvalues of this problem for sufficiently small e —>■ are concentrated on the 
ray 7oo = [—i/\/S,—ioo) and near the segments 7± = [±1, —i/y/S] (see Figure [T]). Of course, 
the density of the eigenvalues is increased as e 0. Based on the method of the paper of 
Morawetz [7j (although problem (|l.ip . (|1.2p was not treated in f^) one can show, that given 
5 > the spectrum of this problem is contained in the (5-neighbourhood of the limit spectral 
graph 

r = 7+U7_, U7oo, 

provided that e < £:o('^)- The author was not acquainted with the paper [2j when he started 
to study this subject. In [13] he described the spectral portrait of problem (jU.l|) . (jU.2|) as 
£ — > based on the properties of fundamental solution of the Airy equation. Moreover, the 
explicit formulae for the eigenvalues near 7_|_ and 7_ were written down. Later on, these results 
were completed and sharpened in the papers of Djachenko and Shkalikov [H], |^. Here we will 
formulate the last results concerning problem ()1.1|) . ()1.2p and sketch basic ideas. 
First, let us note about the following simple fact. 

Lemma 1.1. For any e > the spectrum of problem p.l|) . p.2|) lies in the closure of the 
semi- strip 

n = {A| ImA < 0, -1 < ReA < 1}. 

Proof. Let us consider the operator Ly = iey" + xy generated by boundary condi- 
tions ()1.2j) . Its spectrum coincide with the spectrum of the problem. Obviously, the values 
of the quadratic form 

iLy,y) = -ie{y',y') + {xy,y) 

lie in 11, provided that \\y\\ = 1 and y belongs to the domain S}(L) of the operator. The values 
of the quadratic form form the numerical range of the operator L. Now, the assertion follows 
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-e 



Figure 1. 



from the well-known fact: the spectrum of any operator lies in the closure of its numerical 
range. □ 

The substitution ^ = {—i6)~^^^{x — A) reduces problem ()1.2|) to the following one: 

z{^i) = z{^,) = 0, 

Thus, we have the spectral problem for the Airy equation. This problem is uncommon, since 
the spectral parameter is involved only in the boundary conditions, and the end points change 
depending on the spectral parameter. However, the eigenvalues of this problem can be found 
by standard means. Namely, the eigenvalues are determined by the equation 



A(A) 



^(^i) ^(6) 



0, 



where v{^) and w{^) are linear independent solutions of the Airy equation. Let us take the 
Airy function t>(^) subject the following asymptotic equality 



;i.3) 



^ e An-5, C OO, 



where A^_^ = I I arg^| < tt — 5}, 5 > 0, and the main branches of ^^^^ and ^^^^ are chosen in 
A,r-(5- A domain, where representation ()1.3p is valid, can be enlarged. We take advantage from 
the following result obtained in p,. 

Lemma 1.2. The Airy function f (^) preserves asymptotic representation ()1.3|) in the domain 
A=|e lel>l, |arge|<7r-^|e|-^/^ln|e||, 
i. e. ()1.3|) is valid for all^ G A, and the remainder is majorated by C]^]'^^"^ uniformly for ^ E A. 
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Figure 2. 

Proof. We make use of the well-known representation [U] 



(1.4) v{i) - 

Notice, that in the intersection of A with the second quadrant of the complex plane the estimate 

is valid. Moreover, the first summand in ()1.4|) has representation ()1.3|) in this intersection. 
Using the symmetry principle we obtain the assertion of Lemma. □ 



Let us consider the function 
1 



;l5) 



/(A) 



9 



[(1 - A)3/2 - (-1 - A)3/2] 



which will be used in the formulation of the main theorem of this section. The function /(A) 
is holomorphic in the domain LI where the spectrum is located. We choose the branch by the 
condition /(— z/-\/3) > 0. 
Denote 

A„ = {A G C I I argA| < a} . 

Fix an arbitrary number a > 2~^/^3^^/'' and consider the domain (see Figure E]) bounded 
by the lines ReA = ±1 from the right and from the left, and by the lines passing through the 
points ±1 and the point 



:l6) 



V3 



from the above. The constant 2 ^/■^S ^/'^ for the choice of a is explicit in our calculations, we 
can not make it smaller. 

Lemma 1.3. The function /(A) has the following properties. 
(1) /(A) is holomorphic in the lower half-plane and takes real values on the negative imag- 
inary axis; 
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(2) for all A G the values of f{\) and —if'{\) belong to the sector A^t/q, and the values 
of /"(A) belong to the lower half-plane; 

(3) /(A) increases monotonously and tends to oo as X —too along the negative imaginary 
axis; 

(4) /(A) =2ya + 0(|A|-3/2), /'(A) = VV«A + 0(|A|-3/2) as X^oo, XeD,; 

(5) f{D„)<Ref{X) forXeD^; 

(6) |/(Ai) — /(A2)| ^ """^/^^^ ^ -^o"' where the constant c does not depend on 
Xi, X2; 

(7) the function /(A) takes real values in Dcr only at the imaginary axis. 

Proof. All these assertions follow easily from representation ()1.5p . The details can be 
found in □ 

Let {rfc}^]^ be the increasing sequence of the zeros of the function v{—r), where f (^) is the 
Airy function with representation p.Hj) . Then > 0, and the asymptotic equalities 

3n /, ^ f 1 



(1.7) +0^j^J, *==1,2,... 

are valid. It is not obviuos, that the numeration in these formulas starts from k = 1. The 
explanations are given in [3]. 
Set 

5, = 5^{e) = ae\\Yiel a > 2-'^'^?>-^l\ 

Denote 

yUfc = -1 + e-"'/V/Vfc, fc = l,2,...,fci,fci + l 

/ifc = 1 - e-'''e^/hu, A; = 1, 2, . . . , fci, fci + 1, 

where the number ki = ki{e) is the largest integer such that the inequality e^'^'^r^.^ < 2/v^ — (Jo- 
holds. The numbers /i^ and /i^ lie on the segments 7_ = [1,— i/v^] and 7+ = [—l,—i/y/3] 
respectively. They come up to the 5o--neighbourhood of the knot-point —i/^/3, but only the 
last ones get inside this neighbourhood. 

Next, let us construct a sequence of numbers on the imaginary axis. By virtue of Lemma IT!^ 
f{—ip) is the increasing function for p G M"*". Moreover, 

/K) = /(-^(l/v^ + 5.))>0. 

Given e > choose the smallest integer ko = fco(^) such that f{d^) < irkQE^^'^. Then, each 
equation 

fi-ip) - -nke^/"^ = 0, A; = /cq - 1, A;o, A;o + 1, . . . , 

has the only solution p = Pk > —id^r for k ^ ko and Pko~i < —id^- The sequence {pfej^^^i is 
monotonous and pk oc as k ^ 00. 

Now, we are ready to formulate the main result on the eigenvalue portrait of model prob- 
lem dn}, (113). 

Theorem 1.1. Let a > 2-3/23-3/4 

any fixed number and 6cr = a6\lne\. Set 
^(t) = ^ Re (2e"'/6 - tf'\ te [0, 2/V3]. 
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Denote by the neighbourhoods of the points /i^, 1 ^ k ^ ki, of radius 

7fc = C exp (— (er-'^/^rfc)) , by the neighbourhoods of the points —ipk, k ^ k^ — 1, 
of radius Cp^^e, and by Uq the S^j-neighbourhood of the knot-point — i/-\/3- 

There are numbers C = C{cr) and Eq = Bo{a) not depending on e and k such that for all 
e < So the circles {U^}i^~^^ and {Uj^}^^_i contain the only eigenvalue of problem (jl.lj) . (jl.2|) . 
Moreover, the eigenvalues in are pure imaginary. All the other eigenvalues lie in the circle 
Uo. 

In other words, the pure imaginary eigenvalues in the domain have asymptotics 

(1.9) \k = ~i{pk + epl^Oil)) , k = ko-l,ko,ko + l,..., 
the eigenvalues near the segments 7± have asymptotics 

(1.10) = /i^ + exp(-e-i/V(^'^'^fc))0(l), k = 1,2, ... ,ki,ki + 1, 
and there are 

a2i/233/4|in£| + 0(l) 

eigenvalues inside the circle Uq. In all these formulae the quantity |0(1)| is estimated by a 
constant C depending only on a. 

Proof. Here we outline only the main ideas, the details can be found in 
Let us construct the characteristic determinant of the problem using the Airy solutions f (^) 
and w{^) = V (e"^''*/^^). We have 

A(0 = v{^^)v (6-2-/^6) - V (e-^-^i) ^(6). 

Step 1. First, we investigate the zeros of A(A) in the domain D^j. Notice that the variables 
and e"^'^*/^^^, j = 1,2, lie in the domain A of Lemma fl. 31 provided that A G D^. Therefore, for 
all A G -Do-5 we can use asymptotic representation (jl.3|) for the functions v{^j) and v (e~^'^*/^^j) . 
After simple transformations we get the equation 

which is equivalent in the domain to the equation A (A) = 0. Taking the logarithm of both 
sides, we come to the equations 

(1.11) /(A) -7rfei/2 ^ A-=^/V/2o(l), keZ, XeD^. 

Now, recall the definition of the number ko and the estimate |/'(A)| > C|A|^^/^, A G D^^, which 
follows from Lemma 11.31 Then, using the Rouche theorem, we find that the equation (jl.lll) 
has the only root Xk inside the circle U^, provided that k ^ fco + 1. This root is necessarily 
pure imaginary, since the spectrum of the problem is symmetric with respect to the imaginary 
axis. 

Actually, we can prove the same for the indices fco and kQ — lii we consider a larger domain 
with 2-2/33-1/4 <:(j' <a. It is easily seen due to properties of /(A) proved in Lemma [Ol 

that there are no roots in the domain Dcr outside the circles {U^}^^_i. 

Step 2. While investigating the zeros of A(A) near the segment [—1, —i/\^] it is convenient 

to rewrite the equation A(A) = in the form 
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Let us consider the trapezium fig. C 11, which is bounded by the hues ReA = —1, ReA = 0, 
the real axis and the hne passing through the points 1 and (3^ = -I + 6^^/^(2/73 - 5^). The 
point (3a- is the intersection of the segment 7_ and the circumference of radius 5^ centered at 
the knot-point —i/y/?). The variables ^i, e~^'^*/^^i, e~^'^*/^^2 take the values in the domain 
A of Lemma 11.21 if A G VL„. Hence, we can use asymptotics ()1.3p for the functions involved 

in (jl.l2j) . It is easy to show that the right hand side of (jl.l2p is majorated by C exp ^~|^2^^) • 

Introducing the new variable r = (A + l)e*'^/^ we get from ()1.12|) the relation 

(L13) V{r) := J^^j/.'^!/!,^) = ^(l) exp [s-^'Mr)) , r = (A + l)e-/^ 

The function V{r) has the simple zeros e^/^r^. A; = 1,2, . . ., where have asymptotics ()1.7|) . 
Using the Rouche theorem (here technicalities are omitted) , we obtain that equation ()1.13p has 
simple roots in the neighbourhoods of the points e^^^r^ of radius jk = C exp (— (e^^'^rk)), 
provided that 1 ^ k ^ ki + 1. Coming back to the variable A, we obtain the assertion of the 
theorem about the eigenvalues near the segment 7_. It is easily seen from ()1.13|) that there are 
no other zeros in except the zeros lying in the circles {U^}i^~^^. We pay attention that A^ 
lie in exponentially small neighbourhoods of the points /i^ provided that — i/^/S] > c > 0. 
However, the radii of the circles f/^ increase as /i^ —i/y/3, since (fi{t) ^ as t ^ 2/^3. 
It is easy to calculate that (3k < Ce'^ with some a > 1/2 (depending on a) for all k ^ ki + 1. 

Step 3. We have to show that there are no eigenvalues in the domain H~ \ [Vt^ U -Do-), where 
H~ is the intersection of the semistrip H and the left half-plane. Using representation ()1.3|) it 
is easy to show that for A belonging to this domain the inequality 

Hii)v (e-2-/^6)|>|^(e-^^ni)^(6)| 

holds. 

Step 4- Denote by N{p,e) the number of the eigenvalues of the problem ()1.H) . p.2|) lying 
above the line Im A = —p. Suppose that p > l/VS + 5, S > 0. Fix numbers and e (e < eo) 
and define numbers ko and k as the largest integers satisfying the conditions 

Then 

\N{p,e)-N{p,eo)-{k-ko)\^2. 

Therefore, 

N{p, e) = k + 0(1) = -^f{-^p) + 0(1). 

It is easily seen that N{p,eo) = ko = for sufficiently large eo and 0(1) in the last formula 
takes the values or ±1. 

Now, using asymptotics p.lO|) we can calculate the number of the eigenvalues in the circle 

I 4/2 ^'"^ 



where |0(1)| ^ 3. After simple transformations we obtain 

91/233/4 

Na = 1 In^l + 0(1) + 0(1), 

TT 

where o(l)^Oas£^0 and |0(1)| ^ 3. This completes the proof. □ 
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In addition to this theorem we shall say some words about the movement of the Xkis) as 
e — s> 0. Since the function f{—ip) is monotonous at the semiaxis M"*", it follows from the proof 
of Theorem 11.11 that all the eigenvalues Xkis) lying on the imaginary axis are simple and move 
up as e — >■ 0. Due to the symmetry they can jump off the imaginary axis only in the case when 
Afc gets into the circle Uq and catches up the preceding eigenvalue A^-i. After the collision 
they leave the imaginary axis. We have no information about their behaviour in the circle Uq- 
However, if we proceed to diminish e, the eigenvalues go out from the small circle Uq and move 
along the segments 7_ and 7+ coming to these segments exponentially close. 

2. Model problem: the case of a monotonous analytic function q{x). 

Here we shall show that the spectral portrait of the explicitly solvable model prob- 
lem p.ll) . p.2|) is not accidental; similar phenomena are observed for a wide class of monotonous 
analytic functions q{x). Certainly, another method will be used for the study of general prob- 
lem. Namely, we shall make use of the phase integral method or the so-called WKBJ-method. 
The results of this section were announced by the author |14j . 

We consider the problem 



We pay attention that here for convenience the small parameter is denoted by e"^ instead of e in 
equation (jl.ip . It might be assumed for the simplicity that q{x) admits an analytic continuation 
in the whole complex plane, although only the analyticity in a neighbourhood of the segment 
[— 1, 1] is required in the sequel in addition to the conditions formulated below. 

Let [a,b] be the range of the function q{x) defined for x G [—1,1]. Denote by L{e) the 
operator corresponding to problem ()2.H) . ()2.2|) . Obviously, the values of the quadratic form 
{L{e)y,y) for y G D{L), \\y\\ = 1, lie in the semistrip 



Hence, for any e > the eigenvalues of problem ()1.1|) . ()1.2|) lie in this semistrip (see Lemma IT!^ . 
Let us formulate the main assumptions for q{x). 

(i) The function q{x) is real for x G [—1, 1] and there is a domain G G C such that q{z) 
is analytic in G and maps G bijectively onto H (here the overline implies the closure 
of the domains). 

(ii) For any c G (a, b) the preimage of the ray Vc = {X \ X = c — it, ^ t < 00} 
is a function with respect to the imaginary axis, i. e. any line Im A = const either 
intersects the preimage of the ray only once, or has no intersection points. 

Condition (i) implies that q{x) is strictly monotonous on the segment [—1, 1]. Without loss 
of generality we assume that q{x) is increasing. In this case the domain G belongs entirely 
to the lower half-plane (otherwise, we get a contradiction with condition (ii)). By virtue of 
the symmetry principle, the function q{z) maps bijectively G U G* U (—1, 1) onto the strip 
a < Re A < 6, where G* is symmetrical to G with respect to the real axis. In particular, 
q'{x) > for all a; G (—1, 1). The above conditions are satisfied, for example, for the functions 
q{x) = sin(7ra;/2), g(x) = {x + 1)^ etc. 



(2.1) 
(2.2) 



ie'^y" + q{x)y = Xy, 
y{-l) = y{l) = 0. 



H = {A I Im A < 0, a < Re A < 6} . 
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Figure 3. 



Let us consider the following functions in the closure of the domain 11: 

1 

Q(A) = j ^/iiqix) - A) dx, AeTT, 
-1 

±1 

Q^(A) = ± J V^(g(0-A)c/e, A en, 



where C,\ is the only root of the equation g(^) — A = lying in the domain G. The branches of 
these functions can be chosen arbitrary. To be certain, fix the branches by the conditions 

1 

(2.3) Q(a) = y v/^(g(a?y^^cia; = e^'^/^a, a > 0, 

-1 

(2.4) Q+{a)=Q{a), Q+(A) + Q-(A) = Q(A). 

Certainly, the functions Q, and are analytic in 11 and continuous in 11. 
Define the following curves in the semistrip 11: 

7oo = {AGn| ReQ(A) = 0}, 
7± = {AGn| Reg^(A) = 0}. 

These curves are depicted in FigureElfor the case q{x) = (a; + 1)^/4. Some parts of these curves 
are drawn by dotted lines. We shall show below that these parts play no role in the description 
of the spectral portrait. The remaining parts form the limit spectral graph, and each point of 
this graph accumulates the eigenvalues of the problem as e ^ 0. 

First, we prove some important properties of the curves 7-1- and 700 ■ 



Lemma 2.1. The curve 7+ f7_j passes through the point b (point a) of the semistrip U. 
Both curves are the functions with respect to the real axis. 



SPECTRAL PORTRAITS 



12 



Proof. Consider the curve 7+, the proof for 7_ is similar. Fix a number c G (a, 6). By 
assumption (i) there is the point ^ 1) such that q{C,c) — c = 0. Let \ = c — it, t ^ and 
be the root of the equation g(^) — \ = 0, ^ G- We have 

Q+(A) = y" V^(?(0 - A) rfe = [/+/) 0(g(O-A)rfe=:i^i(A) + F2(A). 

The branch of the function Q'^{X) is defined by condition ()2.4|) . Since Q'''(A) is continuous in 
n, we have 

1 



Notice, that 

1 

ReF^(A) = -- / Re^=^==rfe<0, A = c - zt, 
^ / V (?(0 - c) + 

since < arg ~ c) + zt] < 7r/2. Hence, the function ReF2(c — zt) decreases monotonously, 
moreover, ReF2(c — it) —>■ a.s t 00. Further, denote by ^(/i) the preimage of the segment 
X = c — ifi, ^ fi ^ t, under the map q{z). Due to assumption (ii) we can parameterize this 
curve as follows C,{s) = —is + r(s), ^ s = ^ s{t), ^ 0. 

While calculating the function -Fi(A) we can take the integral along the curve ^(/i). Then, 
Cc 

i'i(A) = j VMO-^)d^= j V^{c-zfL-X)das) = j ,/jr—td{-ts{fi)+r{s{fi))). 

«A s(t) t 

Bearing in mind the choice of the branch, we find that the function 



Re Fi{c — it) = — J y/t — fids{^) 


decreases monotonously as t ^ +00. Therefore, the function ReQ^{c — it) decreases 
monotonously as t changes from to +00. Moreover, this function is positive at zero and 
negative in a neighbourhood of +00. Hence, the equation ReQ^(c — it) = has the only root 
tc < 0, and the curve 7+ is a function with respect to the real axis. Since Q~^{X) — as A — > 6, 
this curve ends in the point b. Lemma is proved. □ 

Lemma 2.2. The curve 700 is a function with respect to the imaginary axis. 

Proof. The branch of Q{X) in H is defined by condition ()2.H|1 . Since a < q{x) < b for 
a; G (— 1, 1), we find that for alH > 

1 



ReQ{b - it) = Re / ^/i{q{x) - b + it) dx < 0. 



-1 
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Further, 

(2.5) ReQ'{X) = -- Re =dx<0, 

2^-1 Vg(a;)-A 

since lm{q{x) — A) > for x G (—1, 1) and A G 11. This imphes that for any fixed t > the 
function Q{c — it) of the variable c G [—1,1] vanishes in the only point. Lemma is proved. □ 

Lemma 2.3. The functions Q~^{X), Q~{X), Q{X) are univalent in the semistrip U. In par- 
ticular, the functions ImQ+(A), lmQ'{X), lmQ{X) are strictly monotonous along the curves 
7+, 7_ and^oo, respectively. 

Proof. It was shown in Lemma I^TD that Re-^Q^{c — it) < for t > and any fixed 
c G (a, 6). Consequently, Re-j^Q~^{X) > for A G 11. This implies the univalence of Q'^iX) 
in n. Since R.eQ~^{X) = along the curve 7+, the derivative of the function ImQ'^(A) along 
this curve does not vanish and preserves the sign. The same assertion is true for the function 
Q~{X). The univalence of the function Q{X) in 11 follows from inequality ()2.5|) . Then, the 
function ImQ(A) is strictly monotonous along the curve 700- Lemma is proved. □ 

Lemma 2.4. IieQ^{X) > if X lies above the curve 7+ and ReQ^(A) < if X lies below 
this curve. Similarly, ReQ~(A) < (> 0) if X lies above (below) the curve 7_. 

Proof. This assertion follows from the proof of Lemma f2. II and representation ()2.4j) . □ 

Lemma 2.5. The curves 7+ and 7_ have the inly intersection point in U. 

Proof. We can view the curves 7_|_ and 7_ as the graphs of negative continuous function 
at the interval (a, b). These functions vanish at the point b and a respectively, hence they have 
at least one intersection point Aq. Suppose, that there is another intersection point Ai. We can 
assume that there are no other intersection points at 7+ and 7_ between Aq and Ai. The parts 
of the curves 7+ and 7_ between Aq and Ai form a closed Jordan curve 7, and the interior of 7 
is a simply connected domain. Consider the function 

<5o(A) = / Vi{q{x) - X)dx, ^0 = 6o> 



«0 

which is analytic in 11. It follows from the definition of 7-1- that the harmonic function R.eQo{X) 
vanishes at both curves 7+ and 7_. In particular, ReQo(A) = for A G 7. By virtue of 
the maximum principle ReQo(A) = in the interior of 7. This implies Qo{X) = const in 11, 
although Q{X) 7^ const. This contradiction ends the proof. □ 

Let Aq be the intersection of 7+ and 7_.. It follows from (j2.4j) that the curve 700 passes 
through the point Aq. By Lemma 12.51 has no other intersection points with 7+ and 7_. 
Denote by 7+, 7_ and 700 the parts of 7+, 7_ and 700 between the knot-point Aq and the points 
b, a and —ioo, respectively. Set 

(2.6) r = 7+U7^U7oo. 

We shall show that F is the limit spectral graph of problem ()2.1|) . ()2.2j) . This implies the fol- 
lowing: any point A G F is an accumulation point for the eigenvalues, while the complementary 
points A G C \ F do not possess this property. Moreover, in the end of this section we shall find 
explicit formulae for the eigenvalue distribution along the curves 7-1- and 700- 
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Figure 4. 

Theorem 2.1. Given r > there exists an Eq > such that all the eigenvalues of prob- 
lem (12.11) . ()2.2j) lie in the r -neighbourhood ofT, provided that e < Eq. 

Proof. We shall use in the sequel the results of asymptotic theory for ordinary differential 
equation. A comprehensive exposition of this theory can be found in books [H]; [I]- 
The function 

z 

s{z,\) = j ^im)-\)di = Q 

ex 

plays an important role in this theory. As before, here ^\ is the root of the equation q{C,) — A = 
from the domain G. By assumption the function q{z) is analytic in a neighbourhood of the 
segment [—1, 1] and admits an analytic continuation in G. Hence, there is a number S > such 
that q{z) is analytic in the domain Q = GU Us[—1, 1], where Us[—1, 1] is the (5-neighbourhood 
of [—1, 1] (see FigurelH). Obviously, the function S{z, A) is analytic with respect to the variable 
A G n and continuous for A G 11, and locally analytic with respect to the variable z & G with 
the branching point ^a- 

Fix a number A = c G (—1, 1) and chose the branch of this function such that Re ^(l, c) > 0. 
For other values A G IT the branch of S{z, A) is chosen by the following condition: the values 
S{l,c) and 5(1, A) are analytically connected. 

For a fixed A G 11 the set 

{^1 ReS{z,\) = 0}, 

determines the lines in the z-plane which are called the Stokes lines. The point ^\ belongs to 
this set and three lines come out from this point. Conventionally, we will use the terms the 
right, the left and the lower Stokes lines. Although these terms are not defined rigorously, the 
identification will be clear from the context. To clarify the situation, we remark that for a fixed 
A G n the domain G may contain some other Stokes lines, besides the Stokes complex with the 
knot-point C,\ (such a line is depicted near the point —1 in Figure |3I). In this connection the 
following assertion plays an important role. 

Lemma 2.6. Let A G 11 \ {7+ U ■j-}, where 'j± are defined in (j2.6j) . Consider all the cases 
(see Figure\^: 

(1) The point A lies above both curves 7+ 7_; 
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(2) A lies under the curve 7_, but above the curve 7+; 

(3) A G 7+ \ 7+, i. e. A belongs to 7+, but lie under the curve 7_; 

(4) A lies under the curve 7+, but above 7_; 

(5) A e 7_ \ 7_, i. e. A belongs to 7_ and lies under 7_|_; 

(6) A lies under both curves 7+ and 7_ . 

Then, in the first case the left and the right Stokes lines intersect the interval (—1, 1) in points 
c~ = c^(A) and = c~^{X), and no other Stokes line intersect the segment [—1, 1]. In the second 
case only the right Stokes line intersect the interval (—1, 1) in a point = c"*"(A) and no other 
Stokes lines intersect the segment [—1,1]. In the third case the right Stokes line intersect the 
point 1 and no other Stokes lines intersect [—1, 1]. The fourth and the fifth cases are similar to 
the second and the third ones, only the roles of the points —1 and 1 are changed. Finally, in 
the sixth case no Stokes lines intersect [—1, 1]. 

Proof. The branch of the function S{z,X) is chosen by the condition ReS'(l,c) > 
for some (and hence for all) c G (—1,1). Viewing S{z,c) as an analytic function on 
z e G Li [-l,c) U (c, 1], we find that ReS{z,c) > for c < 2; ^ 1 and ReS{z,c) < for 
— l^z<c. It is proved in Lemma ITT] that the function ReS'(l, c — it) = ReQ"'"(c — it) de- 
creases monotonously as t runs from to +00 and vanishes in the only point. The Stokes lines 
continuously depend on A, hence, for small values t G [0, to] these lines intersect the interval 
(— 1, 1) in points c~(t), c+(t) which are close to the point c (see Figure IH). For all z G (c"*", 1] 
we have 




Re S{z, c — it) = Re \ / + / A/i(g(x) — A) dx = Re / a/ i{q{x) — A) dx > 



since the point lies on the Stokes line and the values of q{x) — X belong to the first quadrant 
of the complex plane for X = c — it and x ^ c. This inequality shows that no Stokes lines 
intersect the set (c^,l]. Similarly, no Stokes lines intersect the set [— l,c~), and the interval 
(c-,c+). 

Further, the function 

1 

C+ 

decreases monotonously as t grows from 0. Therefore, the point = c^{t) moves to the right 
and reaches the point 1 as X = c — it comes to the curve 7+. Similarly, the point c~{t) moves 
to the left and reaches to the point —1 as X = c — it comes the curve 7_. This analysis makes 
obvious all other assertions of Lemma 12.51 □ 

Now, let us recall an important concept of canonical domain for equation ()2.H) . A domain 
Q\ in the z-plane is called canonical if the function S{z, X) is univalent in fix (we do not define 
here maximal canonical domains). It follows easily from the definition that any domain not 
containing the points of the Stokes graph is canonical. Moreover, Re S{z, X) preserves the sign 
for z belonging to such a domain. This fact implies that any domain containing only one Stokes 
line remains to be canonical. 



Lemma 2.7. Given A G n \ (7+ U 7„) there exist a path connecting the points ±1 and a 
canonical domain Qx which entirely contains this path. 
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Proof. Consider, for instance, the case when the point A hes above both curves 7+ 7_. 
Recall the notation Q = GUUs[—l, 1] and consider the domain where is the domain 

in fl bounded by the left and the right Stokes lines and containing the interval (c~, c+) C (—1, 1) 
(see FigureE}. The domain contains the lower Stokes line outgoing from C,x and, probably, 

some other Stokes lines. However, these other lines do not intersect the lines outgoing from the 
point (this is general property), and do not intersect the sets [— l,c~) and (c+, 1] by virtue 
of Lemma f2. 61 Hence, there is a path in \ connecting the points ±1 which intersects only 
the lower Stokes line, and there exists a neighbourhood of this path containing no Stokes lines 
but a part of the lower one. Such a domain is canonical (as we have noted before). All the 
other cases of disposition of a point A (see Lemma 12. 6|) can be treated analogously. Lemma is 
proved. □ 

The following well-known fact will be essentially used in the sequel. 

Lemma 2.8. Given fixed A G H equation (j2.H) possesses two linear independent solution of 
the form 

(2.7) v^{z, A) = ^ ^ x,ii/4 ^^"''^^''^ (1 + 0±{e)) , 

where the function 0± are subject the estimate 

\0^ie)\^Ce 

with a constant C not depending on z, as z varies on a compact K belonging to a canonical 
domain VL\. Moreover, if K G Qx as X varies on a compact K' in the X-plane, then the estimate 
holds with a constant C not depending on z & K and A G -ft''. 

Proof. See 013]. □ 

Now, we start a direct proof of Theorem 12.11 Fix an arbitrary number r > and denote 
by Ft- the r- neighbourhood of the limit spectral graph F. We have to consider six different 
cases mentioned in Lemma 12.61 for the disposition of the point A with respect to the curves 7+ 
and 7_. Consider, for instance the first case, when A lies under both curves 7+ and 7_. The 
characteristic determinant of problem ()2.H) . ()2.2|) compiled from fundamental solutions ()2.7p 
has the representation 

v-{l,X) v-{-l,\) 
where the function 

r(A) = (%(l)-A)-i/^(^(g(-l)-A))-^/^ 
does not vanish in H. For abriviation we use the Birkhoff notation [1] = 1 + 0{e). It follows 
from representation ()2.8j) that A(A) 7^ if 

(2.9) Re(5(l,A)-5(-l,A)) >0 

and e < eQ = Sq^X). 

Lemma 12.81 asserts that the obtained asymptotic representation for A (A) is valid, if the 
points ±1 can be connected by a path lying in a canonical domain of equation ()2.H) . Lemma ITTI 
guarantees the existence of such a path 7a and a canonical domain ^a. Moreover, the points 
+ 1 and —1 is this canonical domain are separated by a lower Stokes line. Recall that the func- 
tion Re 5(2;, A) changes the sign passing through the Stokes line. We fixed the branch of the 
function S{z, A) by the condition Re S'(l, c) > for c G (—1, 1), and this implies Re 5(1, A) > 



(2.8) A(A) 



T(A) U''~^iS(l^)-S(-l^))y^^^ _ g-e-i{S(l,A)-S(-l,A))^^j 
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if A lies under the curves 7+ and 7_. Therefore, ReS'(— 1,A) < and inequahty ()2.9p holds. 
This inequality remains valid in a neighbourhood Ux of the point A, by continuity. The neigh- 
bourhood Ux can be chosen such small that a path 7 connected the points 1 and —1 does not 
intersect the Stokes lines of the Stokes graph for all ^ & Ux with exception of the lower lines. 
Therefore, the canonical domain Qx containing 7 can be chosen such small that it remains to 
be canonical for all /i G [/a, i. e. we can set = Qx- By virtue of Lemma 12.81 asymptotic 
representation ()2.8|) for A (A) remains to be valid in the whole neighbourhood Ux of the point 
A, and the remainders in the representations [1] = 1 + 0{6) in ()2.8|) can be estimated by Ce 
with a constant C not depending on /i G Ux- 

Fix a number 1 and denote by FI^j the intersection of the semistrip FT with the closed 
circle of radius R with center at the origin. For sufficiently small r the set FIr \ F^ consists 
of three compacts, say, the upper, the left and the right ones. With each point A belonging 
to the upper compact we associate a neighbourhood Ux, that was constructed above, and 
from the cover of by these neighbourhoods take a finite subcover. This procedure allows 
to obtain representation ()2.8|) for all A G K~^, moreover, 

Re(5(l, A) - 5(-l,A)) > c> 0, XeK+, 

with a constant c dependent only on r, and the estimates 10(e) | < Ce hold for the remainders 
in the representations [1] = 1 + 0{e) with a constant C dependent also only on A. Then, we 
obtain A (A) 7^ for all A G if e < and Eq = Eo{t) is sufficiently small. 

Analogously, we can prove the absence of the zeros of A (A) in the left and the right com- 
pacts. It is left to show that the choice of a r-neighbourhood of the curve 700 can be realized 
independently on R, i. e. the eigenvalues with large moduli do not leave the r-neighbourhood 
of the curve 700, vice versa, asymptotically they lie more close to this curve. This fact can be 
proved in the same way as in Lemma 4.4 of paper |18j (see Remark 4.1 in [18^ ). Theorem is 
proved. □ 



Next, we shall obtain an additional information on the eigenvalue behaviour near the curves 
of the limit spectral graph. For this purposes we need the concept of counting function of the 
zeros of holomorphic functions along the curves. 

Let 7 = 7(t), t G [0, 1], be an oriented smooth curve in the complex plane C with end points 
Zq and Zi (the value 00 for these points is admitted). One can order points on the curve as 
follows: Ai -< A2 if Xj = litj) and ti < t^. Denote by 7r(Ai, A2) a curvelinear strip of width 
2r containing 7 as the middle line and having, as the lateral sides, the segments perpendicular 
to 7 at the points Ai and A2. Let a function F{z) be holomorphic in 7^(2:0, -Zi). Fix a point 
Ai G 7 and denote by n{\i, A) the number of zeros of the function F{z) in 7r(Ai, A), if Ai -< A. 
For A ^ Ai define n(Ai,A) = — n(A,Ai). Now, define A^(A) = n(Ai,A) + C, where C is an 
arbitrary constant, as the zero counting function of F{z) in a r-neighbourhood of the curve 7 
(or along the curve 7). Since the eigenvalues of the problem in question are the zeros of the 
entire function A(A), we nay speak in the same context on the eigenvalue counting functions 
along the curves. 
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Figure 5. 

Theorem 2.2. Fix a small number 5 > and denote by ii^ , /x^ and Hk the solutions of the 
equations 

1 

i J V^(g(0-A)de = eirik - 1/4), k e Z, 

-1 



j VMO-^)d^ = - 1/4), kez, 



I 



i 



-1 



Zym^' on the curves 7+, 7_ and 700 (note, that left hand-sides take real values along the cor- 
responding curves). Chose the indices p±, m± and sq such that {/^fc}p!r 
are all solutions of these equations on the curves 7+, 7_ and 700, respectively, lying outside 
the S -neighbourhoods of the points a, b and the knot-point Aq. Then, there exists a number 
C depending only on 5, such that all the eigenvalues lie in the 5 -neighbourhoods Us{a), Us{b) 
and [/^(Ao), and in the circles of the radius Ce^ centered at the points {ti'k}^+-i > {l^k}'p--i 
and {Aifc}^_i- AH these circles contain only one eigenvalue. The eigenvalue counting functions 
along the curves 7+, 7_ and 700 have representations 

iV(A) = — g+(A)+0(l), ^/Ae7+, 
me 

Ar(A) = — g-(A)+0(l), ^/Ae7-, 

N{X) = —Q{X) + 0{1), ^/Ae7oo• 
me 

The remainders in these formulae are estimated by a constant not dependent on A and e, if A 
lies outside the neighbourhoods Us{a), Us{b) and Us{Xo)- 

Proof. Consider, for example, the curve 7+. We shall make use of the transmission formu- 
lae for asymptotic representations of solutions in neighbouring canonical domains. Below we 
formulate the corresponding result in a convenient form, viewing in mind our concrete problem. 

Let A lie in a small r-neighbourhood of the curve 7+ outside some fixed neighbourhoods 
of the end points. Let Cx be the Stokes complex corresponding to the point A, i. e. is the 
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knot-point of C\ (see Figure EJ- One of the Stokes lines, say /i, passes near the point 1 {li 
passes through 1 if A G 7+). As before, denote by VL the union of G and a (5-neighbourhood of 
[— 1, 1]. Consider the canonical domain f2i lying in between the lines I2 and and containing 
the line /i, and the canonical domain VL2 in lying between the lines and li and containing 
12- Let ^^(z. A), vj{z,X) be the pairs of solutions having in canonical domains j = 1,2 
asymptotics ()2.7|) . Asymptotics ()2.7|) for the pair of solutions f^, in the domain fl2, does 
not remain valid in the domain Qi. However, there is a connection formula for the solutions. 

Lemma 2.9. The following transmission formula is valid: for z E fli one has the represen- 
tation 

(2.10) = eW6 f-^l^] [M\ fvti^, A) 



A(A) 



V2 iz,X)J \ 1 J (z, A)^ 

where [1] = l + 0{e) and 10(e) | < Ce with a constant C depending on z and A. However, given 
any compact K in Qi = fli{X) there exists a neighbourhood Ux of the point A such that the last 
estimate holds for all z & K and all A belonging to U\ with a constant C depending only on K 
and U\. 

Proof. See the monographs [S] and [S]. □ 

Now, we can complete the proof of the theorem. Consider the characteristic determinant 

(-1,A) vtHA) 
(-1,A) V2{IA) ■ 

For A lying near the curve 7+ the points —1 and +1 belong to the domains ^2 and fii, re- 
spectively. Using asymptotics (j2.7p and (j2.1(J|) and abreviating by the term e~*'^/^T(A) (see the 
proof of Theorem 12.11) , we find 

[1] exp(£-i5(-l. A)) -i[l\ exp(£-i5(l. A)) + [1] exp(-e-i5(l. A)) 
[1] exp(-£-i^(-l. A)) [1] exp(£-i^(l. A)) 

The branch of the function S{z, A) is defined by the condition ReS'(— 1, A) =: a(A) < 0, if A 
is located in the r-neighbourhood of 7+ outside Usip) and Us{\q). Notice, that Re S'(l, A) ^ as 
r ^ 0. Hence, we can choose a number r > such small that Re 5(1, A) < a(A)/2. In this case 
the term exp£:~^(5'(— 1, A) + ^(l. A)) decays exponentially, while the term exp(— £^-'^S'(— 1, A)) 
grows exponentially as e 0. So, the equation A(A) = is equivalent up to exponentially 
small terms to the equation 

or 

(2.11) cos 1 1 1 v/%(eFA) - ^ J = 0{e). 



A(A) 



A 



Assume that the term 0{e) in this equation equals zero. Then, the roots near the curve 7+ 
are determined explicitly by the equations 

(2.12) -ig+(A) =£(A;7r-7r/4), keZ. 

The function —iQ^{X) is real along the curve 7+ and monotonous (see Lemma l2.3|) . Hence, 
there exist integers p+ and m_|_, such that for all p^ ^ k ^ m_|_ equations ()2.12|) have solutions 
lying on 7+ outside the 5-neighbourhoods of the end points b and Aq. The existence of 
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Figure 6. 

simple roots of perturbed equation ()2.11|) in Ce^-neighbourhoods of the points {/ifc}^"'' can be 
proved by standard means using the Rouche theorem (see details in Theorem 5.1 of |18j ). The 
representation for the counting eigenvalue function N{X) along the curve 7+ can be obtained 
from formulae ()2.12|1 . This can be done in the same way as in Theorem 5.2 of |18j . 

Certainly, the same assertions are valid for the eigenvalues near the curve 7_. The proof of 
the formulae for the eigenvalues near the curve 700 can be realized simplier. Namely, there is 
no need to use the transmission formulae to get asymptotic representation of the characteristic 
determinant for all A lying below both curved 7+ and 7_. For such values of A the segment 
[— 1, 1] lies between the left and the right Stokes lines, and we can use asymptotics ()2.10p 
simultaneously in both points —1 and +1. After simple calculations we obtain that for A lying 
below the curves 7+ and 7_ the equation A (A) = is equivalent to the following one 

sin-Q(A) =0(e). 

le 

Using standard arguments we get the localization formulae for the eigenvalues and the repre- 
sentation for the eigenvalue counting function along 700. Actually, a sharper analysis can be 
carried out. Namely, the formulae 

Afc = /Ufc + eVfc^O(l), = s, s + 1, . . . , 

can be obtained for the eigenvalues near 700 (here |0(1)| is estimated by a constant C not 
dependent on k and e). To get these formulae, one has to use an analogue of Lemma 4.3 from 
paper |18j . Here we omit details. Theorem is proved. □ 



3. The case of Couette— Poiseuille profile 

Profiles of the form g(x) = ax^ + /9a; + 7, a,P,'y G M, correspond to stationary solutions 
of the Navier-Stokes equation, therefore, they have special interest in connection with Orr- 
Sommerfeld problem ()0.1|) . ()0.2|) . Consequently, the study of model problem ()0.4p . ()0.5p for 



SPECTRAL PORTRAITS 



21 



functions q{x) of this form is quite important. Since 



ax +l3x + 'y = a\ x + —— ] +7 — 



2a J ' Aa^' 

the substitution of the spectral and small parameters 

e = v^e', A = v^(A' - 7 + /3V4a2) 

leads to model problem ()U.4p . (jO.Sj) with the function q{x) = {x — P'Y, (3' = (3 /2a. Further we 
write [3 instead of 

The cases = and /3 7^ are different and they have to be treated separately. In the 
first case (i. e. for q{x) = x"^) the eigenfunctions are either even or odd (see for example). 
Consequently, the spectrum of problem ()0.4|1 . ()0.5|) with the function q{x) = x"^ consists of the 
spectra of two problems on the segment [0, 1]: 

, s ie^y" + x'^y" = \y ie'^y" + x'^y" = Xy 

^'^ ' y{0) = y{l)=0, y'{0)=y{l) = 0. 

The first of these problems is solved in §0] the second one can be solved in the same way: the 
replacement of the boundary condition y{0) = by the condition y'{0) = does not change the 
essence of the matter. The method for the solution remains the same. However, the problem 

y'{-i) = y{i) = o 

with non-symmetrical boundary conditions can not be reduced to the study of similar problems 
on the segment [0, 1] with the monotonous function x"^. The solution of ()3.2j) is essentially 
more difficult. The method for the study of this problem is proposed in 18 (formally, the 
Dirichlet boundary conditions are treated in |18j . however, the problem is studied globally on 
[—1, 1], and the proposed method can be used without changes for arbitrary separated boundary 
conditions). It is easy to see analyzing this method, that the form of the limit spectral graph 
r for the function q{x) = x"^ and the formulae for the counting eigenvalue functions along the 
curves do not depend on boundary conditions. Hence, both problems in (|3.1|) and problem (j3.2|) 
have the same limit spectral graph (see Figure IHI) and the same eigenvalue counting functions 
up to multiplication by the coefficient 2. 

Actually, the method of (18l| works to solve the problem 

. . + {x- (3fy = Xy, 

y(-l) = y(l) = 0, 

with 7^ 0. However, the limit spectral graph for this problem takes a more complicated form. 
The same is true for the eigenvalue formulae and the counting eigenvalue functions along the 
curves of the limit graph. 

Let us formulate the main results concerning problem ()3.3|) . Without loss of generality we 
assume, that /3 G (0, 1). Set 

a = (-1-/5)2, 6 = (1-/5)2. 

Repeating simple arguments of Lemma TTT^ we obtain, that the spectrum of problem (j3.3j) lies 
in the semistrip 

H = {A I Im A < 0, < Re A < 6} . 
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Consider the following curves in 11: 

7o = {A G n I arg A = -7r/4} 



7f> 



A G n 



Re j Vz(g(O-A)rfe = 



7a = < A G n 



7_ = <^ A G n 



-1 



Re J V2(g(O-A)rfe = 



-1 



Re 



z(g(O-A)rfe = 



7o 



A G n 



1 



Re y V^(g(e)-A)rfe = 



with g(^) = — It is easily seen that all these curves have no self-intersections (see 
Lemma 12. 5|) . 

Lemma 3.1. The curves 70, 7b and 7_ /iat'e t/ie only intersection point Ai G H. T/ie curves 
7_, 7a anc? 700 intersect in a point A2 G H. There are no other intersection points of all these 
curves. 



Proof. It can be carried out using the ideas from Lemma f2. 51 



□ 



Denote by 70 the part of 70 between and Ai (i. e. 70 = [0, Ai]), by 7b the part of 7b between 
the points b and Ai, by 7_ the part of 7_ between the points A2 and a, and by 700 the part of 
7oo between A2 and —ioo. 
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Theorem 3.1. Given small r > there is a number Eq = EqIt) such that for all s < Eq the 
eigenvalues of problem (|3.3|) lie inside the t -neighbourhood of the set 

(3.4) r = 7oU7aU7_U76U7oo. 

Proof. The first step in proving of tliis theorem is to clarify the geometry of the Stokes 
graphs of the Weber equation (|3.3p . Actually, this work was carried out in [18j . There are tree 
reasons for points A G 11 to be accumulation points for the eigenvalues as e — *• 0. First, at least 
one point +1 or —1 belongs to the Stokes graph Cx of equation ()3.3|) . The set of such points 
form the lines which we call singular. In our case the lines 7a, 7- and 7^ are singular. However, 
not all points of a singular line 7 belong to the limit spectral graph. We have to exclude the 
points /i G 7 possessing the following property: there is a path in the z-plane connecting the 
points —1 and +1 such that it intersect only one Stokes line of all Stokes complexes Cx for all 
A G Us{fi) with sufficiently small 5 > 0. In our case after elimination of such points we get 
the curves 7^, 7_ and %. Second, there are points /i in the A-plane such that the geometry of 
the Stokes graphs is not preserved in any sufficiently small neighbourhood Us{fJ.) (see details 
in |18p . The points of such kind form lines which we call critical. In our case there is the only 
critical fine 70 (for /i G 70 the Stokes graph consists of one complex, while for /i ^ 70 it consists 
of two complexes). Again, we have to exclude the points of critical lines which possess the 
property that we described before. Then, in our case we get the curve 70. Finally, the curve 



7o 



Re / a/ i{q{x) — A) dx = 



-1 



we call the main line. The part of this line between — ioo and the first intersection point with 
singular or critical lines has to be included in the limit spectral graph. 

This introduction gives an understanding, how to prove Theorem 13.11 Here a complete 
proof is omitted. The details can be found in |2Qj. □ 

Theorem 3.2. The set fl3.4j) is the limit spectral graph of problem fl3.3p . i. e. the points 
A G r and only these points are the accumulation points of the eigenvalues as E 0. Given 
6 > there are numbers Eq = Eo{S) and C = C{5) such that for all e < Eq the eigenvalues of 
problem fl3.3|l lie in the set 

Us = Us{Q) U Us{a) U Us{h) U Us{\i) U [/^(Aa) 

and in the C e"^ -neighbourhoods of the points iik which are determined by the equations 



^ I A/^(g(0 - ^) dC. = e{Txk - 7r/4), /i^ G 75, ke Z, 

-1 

^ / A/«(g(0 - A) rf^ = E^nk - 7r/4), ^ik^la, k E Z, 

-1 

i I 0(g(O - A)c?e = e(vrfc-7r/4), /ifc G 7-, k e Z, 
1 

« / a/^(?(0 - ^)d^ = ^tt/c, Aifc e 7oo, kez, 
-1 

A^ = {2k + \)Ee-'^l\ /ifc G 70, k^ Z. 
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Figure 8. 

The C e'^ -neighbourhoods of all points G F \ f/5 contain only one simple eigenvalue. The 
counting eigenvalue functions along the curves of the graph F have representations 

1 

NiX) = —^ [ V^(g(e)-A)rfe + 0(1), AG 76, 



tire 

-1 



N{X) = — [ ^/^{q{0-X)d^ + O{l), A G 7a, 



-1 



N{X) = — [ v/^(^(ePA)ci^ + 0(l), AG7- 
ine J 



1 



N{X) = — [ VKq{0-X)d^ + 0(1), A G 7o 

tTTE J 



-1 



iV(A) = ^e^^/'A, A G 70- 



Proof. The details of the proof of this theorem can be found in |20j . The spectral portrait 
of problem ()3.3|) for the function q{x) = {7x — 1)^/64, = 5000 is depicted in Figure [7| □ 



4. The Orr— Sommerfeld problem 

The description of global behaviour of the spectrum of the Orr-Sommerfeld problem as 
i? — > cxD have been carried out only for the functions q{x) = x (Couette profile) and q{x) = x^ 
or q{x) = 1 — x^ (Poiseuille profile). Here we point out the papers [1], [3], [7] and |19j . 
Actually, the methods developed in |19j and |20j can be used without essential changes to 
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Figure 9. 

treat the Couette-Poiseuille profiles q{x) = ax^ + /5x + 7. Here we shall only formulate the 
result obtained in [3] and the generalization of the result |19j : the detailed proof requires a 
serious work. Figure |H1 (it is borrowed from the paper f^) gives an illustration of the first 
theorem of this section. Here the spectrum is calculated for a = 1 and R = 4000. 

It was shown in Introduction, that in the case of the Couette profile q{x) = x the Orr- 
Sommerfeld problem takes the form 

—iez" + {iea'^ — x)z = Xz, e = 1/aR, 

^^'^"^ y sinh[a(l - t)] z{t) dt = j cosh[a(l - t)] z{t) dt = 0. 

-1 -1 

To describe the eigenvalue behaviour of this problem as £ — > 0, we consider the rectangular 
coordinate system {t, 7} in the A-complex plane taking the point —1 as the origin and directing 
the axis t along the segment [—1, (see Figure Oj). 

Set 

^Isinh (a (2 -e-*^/^t))| 1 / / ■ u ^^ 

cit) = 2v^l ^, m = ^ argsinh {a (2 - e^/^t)) , 

where the main branch of the argument is chosen. 
Consider in the {t, 7} plane the curves 



e'/\ c 



(t)t3/4 



7±W = ±^ln^^, t>0, 
and fix at these curves the points 

where 
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Figure 10. 

and the indices k^, ki are chosen in such a way that 

£^/^|lne| ^ tfc ^ 2/V3- - f - j £^/2|lne|, k^^ki^ki. 

Theorem 4.1. Denote by the neighbourhoods of the points jji^ of radius 6^ = Ce^^'^t^^^^ 
and by the symmetrical reflections ofU^ with respect to the imaginary axis. Denote by U±i 
the e^^^l \n e \ -neighbourhoods of the points ±1 and by Uq the 2/3 ■ (3/4)^/^6"'^/^ | \n el-neighbour- 
hood of the knot-point —i/\/?). Then, there are numbers C > and £o > such that all the 
eigenvalues of problem 1)4.11) located near the segments [±1,— i/v^] He- inside the circles U±i, 
Uq and {U^}^^^\, {Uk}%X^i! provided that e ^ Sq. The circles and contain only one 
simply eigenvalue. All the other eigenvalues lie at the imaginary axis below the point —i/^/S 
and have representation 

Xk = -i{pk + £0{1)), k = ko,ko + l,. . . , 

where the numbers pk and ko are the same as in Theorem \l.l\ 

Proof. As in Theorem II. H special properties of the Airy functions are explored in the 
proof. Details can be found in □ 

Theorem 4.2. Given r > there is e = eo{T), such that for e < sq all the eigenvalues of 
the Orr-Sommerfeld problem (jO.ll) . ()0.2|) with the Couette-Poiseuille profile q{x) = {x — (3Y , 
(3 G (—1, 1), lie in the r -neighbourhood of the limit spectral graph T of the corresponding model 
problem. The main terms of the counting eigenvalue functions along the curves of the graph T 
have the representations given in Theorem \3.1\ 

Proof. In the case /3 = the proof is obtained in |20j . For /5 7^ the proof remains 
essentially the same, provided that the analysis of the model problem with q{x) = [x — jS)"^ is 
carried out (see [20j). The proof of the results about the counting eigenvalue functions uses the 
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tauberian technique developed in |14j . Figure ITUl shows the spectrum of the Orr-Sommerfeld 
problem for q{x) = x"^, a = 1 and R = 3000. □ 
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